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Abstract 

We study the limit of quasilocal energy defined in [7] and [5] for a fam- 
ily of spacelike 2-surfaces approaching null infinity of an asymptotically flat 
spacetime. It is shown that Lorentzian symmetry is recovered and an energy- 
momentum 4- vector is obtained. In particular, the result is consistent with the 
Bondi-Sachs energy-momentum at a retarded time. The quasilocal mass in [7] 
and [8] is defined by minimizing quasilocal energy among admissible isometric 
embeddings and observers. The solvability of the Euler-Lagrange equation for 
this variational problem is also discussed in both the asymptotically flat and 
asymptotically null cases. Assuming analyticity, the equation can be solved 
and the solution is locally minimizing in all orders. In particular, this produces 
an optimal reference hypersurface in the Minkowski space for the spatial or 
null exterior region of an asymptotically flat spacetime. 

1 Introduction 

This is a continuation of |9] in which the spatial limit of the new quasilocal energy 
defined in |7] and [8] is analyzed. In the present article, we address the question of 
the null limit in Bondi-Sachs coordinates for an asymptotically flat spacetime. Let 
be a spacetime with metric gap in Bondi-Sachs coordinates given by 

-UVdw'^ - 2Udwdr + aabidx" + W''dw){dx^ + W^dw) a, 6 = 2, 3 

where 

W" = 0(r-2), 
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Y'i _|_ V2 

V =\ + o[r-^) 

r 



and the metric (Tab is given by 



^2 + 2Xr + 2(^2 + y2) _2Yr sin 6* 

-2Frsin^ sin^ - 2Xr + 2(^2 + F^)] 



with 



det (Tab = ''^^ sin^ 6*. 
The inverse of the metric g^p is 

^ gwa ^ Q^^«,r _ = [/" V, = [/"^H^", and = CT'^^ 

Throughout the paper, coordinates are labeled by = w,x^ = r, = 6,x^ = (p 
and the indexes are for a, 13, j ■■■ = 0,1,2,3, i, j, k ■■■ = 1, 2, 3, and a,b ■■■ = 2,3. 

At a retarded time w = c, the Bondi-Sachs energy- momentum vector ([1] [6]) is 
defined as 

(S,Pi,P2, A) = — ( /" 2mdS'^, [ 2mXidS^, [ 2mX2dS\ [ 2mX^dS'^) (1.1) 

where m = m(c, 9, </>) is the mass aspect function in the expansion of V, and Xi, 
i = 1,2,3 are the three eigenfunctions sin 6 sin (j), sin 9 cos (j) and cos 9 of the Laplace 
operator A on S'^ with eigenvalue —2. 

We recall that given a spacelike 2-surface S in a spacetime, a quasilocal energy 
i?(S, X, To) is defined in [7], [8] with respect to an isometric embedding X : S — )■ M^'^ 
and a constant future timelike vector Tq G M^'^. For a family of surfaces and a 
family of isometric embeddings X,. of S,, into M^'^, the limit of £'(Sr,Xr,To) is 
evaluated in [HI Theorem 2.1] under the assumption that 

liin = 1 (1.2) 



r— >oo 



where if and Hq are spacelike mean curvature vectors of S,, in and the image of 
Xr in M^'^, respectively. In fact, the limit of E{T,r, X^., Tq) with respect to a constant 
future timelike vector Tq G M^'^ is given by 



lim — 

r-s>oo Svr 



-{%. ^)(iff„i - \H\) - (v-r ^. ^> + (v^,^, A) 



(1.3) 
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where r = —{TQ,Xr) is the time function with respect to Tq, and Jo and J are the 
future timehke normal vectors dual to Hq and H. This expression is linear in Tq and 
defines an energy-momentum 4-vector at infinity. 

In this article, we consider a family of 2-surfaces on a null cone w = c as 
r goes to infinity in Bondi-Sachs coordinates. The limit of the quasilocal energy is 
first computed with respect to isometric embeddings Xr into which are essentially 
unique and satisfy (11. 2p . We show in particular, 

hm^ / i\Ho\-\H\)dK = E, and lim ^ / (V^^^, -^)ciS, = (1.4) 

where [Xi, X2, X^) are the coordinate functions of the isometric embedding Xr into 
R^. We remark that exactly the same limit expression on coordinate spheres of 
asjTiiptotically flat hypersurface gives the ADM energy-momentum in [S] . The com- 
putation is stable with respect to any 0(1) perturbation of Xr in M^'^ and is equiv- 
ariant with respect to Lorentzian transformation acting on Xr- 

In [7] and [8] , the quasilocal mass of a 2-surface E is defined to be the minimum 
of E{T,, X,To) among all admissible pairs (X, Tq) and the Euler-Lagrange equation 
is derived for an optimal isometric embedding. In the last section, we show that 
an analytic solution of the optimal isometric embedding equation can be obtained 
as an 0(1) perturbation of embeddings into a boosted totally geodesic slice in R'^'^ 
whose timelike normal is in the direction of the total energy-momentum 4-vector. 
This solution locally minimizes the quasilocal energy. 

Brown-Lau-York [2] and Lau [1| compute the null limit of the Brown- York energy 
and we compare our calculation with theirs in the following: 

1) Brown- York mass is gauge dependent. After fixing a reference isometric em- 
bedding (either to flat R^ [2] or to the null cone in R^'^[l] ), a gauge is chosen 
arbitrarily so that the limit of the mass coincide with the Bondi mass. In contrast, 
in our case, once a reference isometric embedding is picked, the quasilocal energy 
is determined by the canonical gauge condition (Eq (1.1) in [9]). Our calculation is 
robust with respect to the choice of reference isometric embedding. In particular, 
the reference family can be arranged to be asymptotically fiat or asymptotically null 
in R3'1. 

2) In [2], the momentum part came from the smear energy while in our case, 
the momentum part came from the connection one-form associated with the mean 
curvature gauge. This one form gives the right momentum contribution in the asymp- 
totically fiat well (see [9]). 

3) In [2], the energy and momentum are defined separately. In our case, the 
Lorentzian symmetric is recovered at infinity and the energy-momentum form a co- 
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variant 4-(co)vector. We show that this (co)vector is equivariant with respect to the 
reference isometric embeddings into M^'^. 

Acknowledgement: Part of the work is done while the authors are visiting the 
Taida Institute for Mathematical Sciences in Taipei, Taiwan. 



2 The geometry of 2-surface in Bondi coordi- 
nates 

Let be an asymptotically fiat spacetime with Bondi-Sachs coordinates. Let 
be the 2-surface defined hj w = c and a fixed r. In this section, we compute the 
mean curvature vector H of S,, in and the connection one-form of the normal 
bundle of in the mean curvature gauge. Denote Wa = CTabW' and let 6°'Wa be the 
divergence of the 1-form Wa on Ur with respect to the induced metric aab- 

Lemma 1. Let be the 2-surface defined by w = c and a fixed r. The mean 
curvature vector H of in N is given by 



U 



r dw dx°- r dr 



~ abdrc 

T^w r 
abdwrQ 



J- ab9rd(^ - J- abQwdQ 



(2.1) 



In particular, H is spacelike when r is large enough with 

\H\^ = -^{- + 6''Wa). (2.2) 
Ur r 

Suppose J is the future timelike normal vector dual to H, then 

(V a J, H) = + S'^Wa) - 4^(- + 6''Wa)cT,kdrW''. (2.3) 

rU r rU r 

Proof. By definition, we have 

H = a^\V^4-b-i^^irjf) 

= ^'''ra^Yr + ^^'d^ + ^^^^ ~ ^^^^d^^' d^>'^^d^^^- 

The last coefficient can be computed explicitly as 

pc /Y7 ^ ^ \„dc 

^ab (V_|^^,^)a 
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Thus 

The relevant Christoffel symbols of gai3 are given by 
and 

Kb = -lu-'[d,Wa + daW, - d^aab + Vdra^b " 27,^,1^,] 

where is the Christoffel symbol of the metric aab- When tracing with cr"^^, we 
notice that 

(y^^dc^aab = da (In det aab) = da ln(r^ sin^ 6) . (2.4) 

Thus, we obtain equations (12 .ip and (I2.2p . To compute the connection one- form, we 
rewrite equation (12. ip as 



r dr r dw dx'^ 

V 3 2 3 3 V 3 

= — (— + )— + -(— — — W^—— — — — ) 
r dr r 3w dx^ 2 3r 

where and — W^j^r — are null vectors. 

or aw ox'^ 1 or 

Thus we have 

V 3 2 3 3 V 3 

^ r ~^ 3r r^3w 3x'^ 2 3r^ 

For simplicity, let's denote ^ and ^ — W^-^ ~ and the coefficients 

^ + S^'Wa and ^ by x and y in the following computation. Then, 

(V 8 J, H) = [/^^(V a xfii + yn2, xui — yn2) 

= U''^{[{3bx){-y) + {3by){x)]{ni,n2) - xy{{Wani,n2) - (V a ^2,^1))} 



U [{3bx){-y){ni,n2) - 2xy{V _a_ni,n2) + xy3b{n2,ni)] 



(2.5) 



On the other hand. 



3 3 3 V 3 



because = by direct computation, and is perpendicular to the null normal 
& - W^'^afd - -jf • Substitute in Tl = ^U-^hu - ^U~^abcdrW^, and we obtain 

(V^ni, na) = -^dbU + ^abcdrW. (2.6) 

Plug (12 .op and {ni,n2) = —U into (I2.5p . we derive 

(V^ J, =[/"22/[(<9bx)f/ - xabcdrWci 

and (I2.3p follows in view of the definitions of x and y. □ 



3 Limit of quasilocal energy 

In this section, we compute the limit of quasilocal energy with respect to a family of 
isometric embeddings Xr of as an 0(1) perturbations of a boosted totally geodesic 
slice in M^'^. First we quote the following lemma whose proof can be found in [3]: 

Lemma 2. Let cr^^ be a family of metrics on — S'^ with cr^^ = r'^dab + 0{r) in 
which dab is the standard round metric on S'^. Let Xr = {Xi, X2, X^) be the isometric 
embedding into for r large and Hq be the mean curvature of Xr ■ Then 

\Ho\ = - + 0{r-') and [ |//oMS, = Anr + ^''''^^-^ + 0{r~'). 

We note that up to an isometry of M'^, Xr can be arranged so that the coordinate 
functions satisfy X^ = rXi + 0(1). 

Theorem 1. Let S,. be the 2-surface defined by w = c and a fixed r in an asymptoti- 
cally flat spacetime with Bondi-Sachs coordinates. Suppose Xr is the ( unique ) family 
of isometric embeddings of S,. into M? for r large, the limit of quasilocal energy with 
respect to Tq = (^1 + |aP, a\ a^, a^) is 

lim ^(S,,X,,ro) = — / 2m{^/l + \a\^ + a'Xi)dS^. (3.1) 
r^oo Sit J g2 

Proof. Let {0,Xi,X2,Xs) be the isometric embedding Xr of into C M^'^. In 
this case, is simply the vector (1, 0, 0, 0). By the assumption on aab we can apply 
Lemma [2] and 

/ \Ho\dEr = 87rr + 0{r'^). 
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On the other hand, from equation ( 12.2^ and the expansion for V, we obtain 

2 9m 

\H\ = -- — + S'^Wa + 0(r-3) (3.2) 

and thus 

I \H\dJ:r = 87ir- I 2mdS^ + 0{r^^). 
Next we compute the physical hamiltonian 

V Z-lf I I I I Z-lf I I I I 

From equation (12. 3p and the asymptotic expansions of V and W"" , we derive 

9 9m 4 

(V^ J, H) = -[d,{5-Wa - -^)] + -M, + 0(r"^). 

Let V denote the connection one-from (V^ j^, j^). From the above computation, 

d^v^^V = i[i(A + 2WW. - + + 0{r-% 
The hmit of (SvXijk\' ]^)^^r- as r — > 00 is thus the same as 

hm / Xidw^^Vd'Lr= I Xi\^{^ + 2){5^Wa-2m) + 2m]dS'^ = I Xi2mdS^ 
'■^"oJe, Js^ 2 Js2 

In this case, the reference Hamiltonian term is zero as 7^ is a constant vector. 
In view of expression (II. Sp . the theorem is proved. 

□ 

Next we show that the limit of the quasilocal energy is invariant under any 
0(1) perturbations of embeddings into totally geodesic and that it is Lorentzian 
equivariant. 

Corollary 1. Suppose Xj. = {tq, Xi, X2, X3) is a family of isometric embeddings of 
Tir into M^'^ with tq = tq + 0{r~^) for some function tq on S"^ . Then we still have 

lim Ei^Lr, X„ To) = / 2m{^/TTW + a'Xi)dS^. 
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Proof. Let Xr be the embedding of by projecting Xr onto M.^ which is given by, 
{0, Xi, X2, X3). It is not hard to check that the induced metric by the embedding 
Xr agrees with the standard round metric of radius r up to the top order term and 
its area agrees with that of the standard round metric of radius r up to the second 
order term. The mean curvature of the embedding Xr is then (0, AXi, AX2, AX3). 
By Lemma |2l the mean curvature Hq satisfies 

\Ho\ = '^ + 0{r-^), and j \Ho\dLr = ^nr + 0{r-^). (3.3) 

The mean curvature Hq of Xr is given by 

(Aro,AXi,AX2,AX3). 

The difference between aah and aah is of order 

cTab - = Oil) and - a'^^ = ©(r"^). 

As a result, the difference between the two laplace operators is of order 

KXi- AXi = 0{r-''). 

Hence |^oP - l^oP = 0{r~^), and thus |^o| - l^ol = Oir^"^). 

By equation fll.3p . the hmit of quasilocal energy with respect to the embedding 
{rQ,Xi,X2,X^) is thus 

^ hm / {\Ho\ - \H\)dK = ^ lim / (|i^o| - \H\)dK. 

Unhke the previous case, is no longer a constant vector for such an isometric 
embedding Xr. However, the asymptotic expansion 

l-"o| 

is valid and the energy component is the same as the limit of quasilocal energy of 
the isometric embedding into in view of (13. 3p . 

Next we compute the physical hamiltonian. Since the induced metric on the 
projection still agrees with the standard one up to lower order term, up to an isometry 
of M'^, = rXi + 0(1). The corresponding time function r is 



T = a^Xi) + ro /l + 5^ a2 = -(^ a,X,)r + 0(1). 

i V * * 

Thus the physical hamiltonian remains the same. 

Lastly, we claim that the reference hamiltonian goes to as r goes to infinity. 
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Lemma 3. Let Cr = r^o" + 0(r) be a family of metrics on S^. Given an 0(1) time 
function tq, let be the images of the isometric embedding Xr into M^'^ determined 
by Tq. Let Vq be the vector dual to the one form (V(.)|^, ^jy) on S° then 

div^^V^ = ^A(A + 2)ro + 0(r-^). 

Proof. We need to compute up to 0(r~^). For this purpose, it is enough to 

assume that the embedding is (tq, rXi, rX2, rXs). Using ^^dcXid^Xi = a^c, we 
derive that a normal vector is 



(1, -daTodbXid- 



ah\ 



where a is the standard metric on 5"^. As XidaXi = 0, we check that 

r \Llo\ 

is a normal vector perpendicular to Hq. Thus, which is the unit normal per- 
pendicular to Hq is, up to lower order, given by the same expression. As a result, we 
compute 



i 

The right hand side equals to 

J2kdbro){dad,Xi)a''AX, 



-{d,Tod,X,)a''^ + ^AX, 



\Ho\ 



daAro + 0(r- 



AX,. 



Here one uses again that Yli ^ida^i = and thus from the first term, one has non- 
zero contribution only when the derivative da falls on dcXi- For the second term, the 
leading term of \Hq\ is independent of 9 and (j) and ^^(AXj)^ = \Hq\'^ up to lower 
order terms. Thus one only has contribution when the derivative hits Atq. Direct 
computation using J2i XidadcXi = —cFac shows that 



Y,-i9kTQ)idadcX,Ya''AX, = 



As a result, 



div^^Vo = 77^A(A + 2)ro + 0(r' 



2r 



□ 
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Using the above lemma, the reference hamiltonian at infinity is 

" 1 
'S2 2 



hm / Xidivj^^Vod^ = [ lXiA{A + 2)TodS^ = 



□ 



Corollary 2. Suppose X'^ is another family of isometric embeddings of ar into 
M.^'^ such that X'^ = L^X^ for some X^ in the previous Corollary, and a family of 
Lorentzian transformation Lj. such that the limit of the 5*0(3, 1) part of Lj. converges 
to an Loo, then the energy-momentum ^-vector also transform by Loo- 

Proof. Both \Ho\ and the connection one form j^) invariant under 

Lorentzian transformation, while (Tq, and Vr = — V(To,Xr) are Lorentzian 
equivariant. □ 

For example, if we take a family of isometric embedding X^ into and define 
X'^ = Xr + r, it is not hard to see that the hypersurface spanned by X', is asymptot- 
ically null. 



4 Optimal embedding equation 

The optimal embedding equation for minimizing the quasilocal energy is derived in 
[H Proposition 6.2]. The equation reads 

-iHa''^-a'''a^'^h,d) J^^J^ + divj^i , coshglgl - V6'- V) = (4.1) 

where sinh^ = T^"^ and V is the connection one-form (V^m-dt, rm)- To solve 

I-h'Ia/i+|Vt|2 ^ (■) l^r 1^1' 

for this equation, we start with data on the 2-surface S given by [pah^ y\ Take 

a function r on S and consider the isometric embedding X : (S, a) — M.^ with 

the metric dab = (^ab + TaTb. H and hab are the mean curvature and the second 

fundamental form of S, the image of X in M^, respectively. 

Let X : E — )■ M^'^ be the embedding of the graph of r over £ in M'^'^, and Sq be the 

image of X with the induced metric isometric to a. The optimal isometric embedding 

equation can be written in terms of the geometry of Sq. In fact, the quasilocal energy 

of So with respect to itself as a reference is zero and thus minimizing. Equation (14. ip 

is automatically true on Eg and we deduce 

- {Ha'^' - a--a»%,)^^^ + div^j , ^] = cosh^ol^^ol - V^o - V^,) = 

^yl + |Vr|^ ^yl + |Vt|^ 

(4.2) 
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where Vq, Hq, 9q are the corresponding data on Sq. Equation (14.2p can be checked 
directly for spacehke surfaces in R'^'^. 

Subtracting equation (14. ip from equation (14.21) . equation (14. ip is then equivalent 

to 

div^[ '^^ =(coshg|ij"| -cosh6>o|/7o|) - V(g-go) -V + Vo] = 0. (4.3) 
a/1 + |VrP 

By the definition of 9 and 6*0, we derive 



cosh -cosh ^ol^ol = W|^P + -^^^^ - Wl^^or ' ^^^^^ 



l + lVrP \ ' l + |Vr| 



and 



|//||iJo|\/l + iVrP'V ' ' 1 + |Vt|2 y l + |Vr| 



Set 



/ 



1 + |Vt|2 







2 _ 


l^ol 


2 


V\Hr\ 


P(l + 


|Vr| 


P) + (Ar) 


+ 


\/l^oP(l + 


|Vr| 


|2) + (Ar)2 



(4.4) 



and equation (14.10 is equivalent to 

div^ifVr) - A[sinh-^( ,f'[£ , )] - (dzi-sT/ - dzt;^^)) = 0. (4.5) 

l-n ll^ol 

In this equation, V are \H\ come from the physical data, and Vq and \Ho\ only 
depend on the embedding X : E — >■ M^'^, while r = — To) depends on both X and 
Tq. Equation (14. 5 p together with the isometric embedding equation {dX, dX) = a 
form the optimal isometric embedding system. 

We shall solve the system for a family of spacelike 2-surfaces at null or spatial 
infinity such that the family of isometric embeddings X^ into M^'^ is of the form 
Xr = BrXr where is a family in 150(3, 1) and Xr is an 0(1) perturbation of 
isometric embeddings into M^. 
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We observe that momentum become an obstruction to solving the optimal em- 
bedding equation for Xr and then discuss how this can be resolved by boosting the 
embedding by Br- The discussion covers the spatial infinity case discussed in [9] 
as well. In the last subsection, we show the solution obtained is locally energy- 
minimizing up to lower order terms in r. 

4.1 Embedding near 

In this subsection, we study the geometry of a family of isometric embeddings 
that is near a totally geodesic in M'^'^. 

Lemma 4. Suppose Xr is a family of isometric embeddings into M.^'^ for a given 
family of metrics = r'^a + ra^^'^ + YlT=o r~''a^~''^ and 

oo 

Xr = rX(i) + r-^X^~^^ (4.6) 

fc=0 

with X^^^ = X = {0, Xi, X2, X^) , the standard embedding of into M^. Denote by 
f^^^ the time component of X^'^'^ and by G{f^^\ ■ ■ ■ ,f^~^\a) a term that depends on 
f(°),-- - and a. Then 

00 

k=3 

where Hq = G{f^^\ ■ ■ ■ , f^~'^~^^\ a) for k > 3 and depends only on a. 

Proof. In the proof, we suppress the subscript r and write X for Xr and a for ar- 
From the expansion oi a, a = r'^a -f rcr'^^^ -|- ■ ■ ■ , we deduce 

i=l 

This can be transformed into a linear elliptic equation for xf^ which can be solved 
(see section 6 of [S]). In general, 

i=l m=0 

and thus x\~^~^'^ is determined by f'^°\ ■ ■ ■ , f^"') and a. 
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Recall the mean curvature is given by Hq = AX where A is the Laplace operator 
with respect to a. For a function g on S"^, we compute 

Ag = a'^\dad,g-j:,d,g) 

= (r-^a'^^ + r-V(-3)'^^ + 0{r-^)){Vl,g - r-Wab^^"d,g + 0{r-^)) 

where = 7^^ + r~^^ab^^^ + ■ ■ ■ is the expansion of the Christoffel symbol 7^^ of 
(Tab and Vli^g = dadbg — Jab^cg is the Hessian of g with respect to aab- Therefore, we 
obtain the following formula: 

Ag = r~'Ag + (a^^'^^'Vlg - a'^^i^'^' d^g) + 0{r-') (4.7) 

By (14. 6p . Hq = AX has the following expansion: 



Ho = rAX(^) + ^r-'=AX( 



fc=0 

Since X^'^'' = X, the standard embedding of S"^, we compute 



where we use V^^X = — cra^X. 
Therefore, we obtain 



00 

Ho = AX = -2r-'X + r~^Ht'^ + ^ r-'nt'^ 

k=3 

where 

= AX(°) - (aa,a(-^)'^^X + a'^S(-^):,9.X) 

and 

Hf^ = AX("^-+2) + ^(f^o), ■ ■ ■ f(-'=+3)^ ^) 

for k > 3. 

We compute 

00 3 

|Fo|^=4r-2+r-3(2AX,AXf^+G'(a))+^r-'=-i(2^AX,AX^''+')+G'(f^^^ , f(-'=+=^), a)). 

fc=3 i=l 

Since X-*-"'^'*'^'' depends only on f^''^ ■ ' ' , f^"'^"'"'^) and X^^'^ depends only on a, we 
obtain 

00 

\Ho\^= 4r-2 + r-^G{a) + ^ r~''~'G{f^''\ ■ ■ ■ , f^-'^), a) 

k=0 

and the expansion for \Hq\ follows by taking square root. □ 
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Lemma 5. Under the same assumption as Lemma \^ Let Vq be the one form 
(^(■)J7^' ^15]'^ ^"^ag'e of Xr then 



r V, 







^ oo 

divj:^Vo = r-3[-A(A + 2)f(°)] + 

k=A 

where 

v^-^) = iA(A + 2)f(-^^+=^) + G(f ■ ■ ■ , f(-^+^), a). 

Proof. To compute the connection one-form in mean curvature gauge, we need a 
normal vector I that is perpendicular to Hq. Suppose / is of the form 

oo 
1=1 

with J(°) = (1,0,0,0). The condition {I,daX) = implies 

(j(°),a,xW) + (/(-i),a,l) = o, 

and (/, Hq) = implies 

{I^'\Ht'^) + {I^-'\-2X) = 0. 

Since X, diX, d2X form an orthonormal basis in M'^ and we can assume is 
perpendicular to it is not hard to check that 

j(-i) = a,f(°)a^^a5l - -Af(°)x 

and in general 

j(-fe) = a^f - ^Af^-'^+i)^ + G(f(°), ■ ■ ■ , f(-'^+2)) 

for A; > 1. 

The connection in mean curvature gauge is thus 

^-{1,1) \J^o\ ^-{I,I)\Ho\ 
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We compute 

oo 

{VaI,Ho) = -2r-\daI^~'\X) + J2r~\-2{daI^-''-'\X) + ■■■ + {daI^-'\ Ht'^'^)) 

k=3 

oo 
fc=3 

Since the leading term of \Hq\ is |, the leading term of (/, /) is —1, and the 

leading term of div^^a is r~^diva for a one- form a on S"^, we obtain the desired 
expansion for div-^^Vo. 

□ 



4.2 Boost in ^ 

Suppose is a family of spacelike 2-surfaces in spacetime such that 

(1) The induced metric satisfies ar = r^a + 0(r). 

(2) The norm of the mean curvature vector satisfies \H\ = - + + 0{r~^). 

(3) The connection one-form in mean curvature gauge V satisfies div-£^V = — h 
0(r-4). 

These assumptions hold on coordinate spheres of an asymptotically flat hyper- 
surface as well as the r level surfaces at a retarded time in Bondi-Sachs coordinates. 
Altogether they guarantee the limit of the quasilocal energy- momentum (e,pi,p2iP3) 
with respect to isometric embeddings of cr^^ into is well-defined. 

We assume the family of isometric embeddings Xr : T, x [ro, oo) — )■ M^'^ is of the 
form 

Xr = BrXr (4.8) 

where X,. is a family of isometric embeddings of into M.^'^ that is analytic in r 
and i? is a family of elements in iS'0(3, 1) that is analytic in r. We assume that 
B = linifc^oo Bk where 

B, = eT^''-''---e'r^'~\'''''\k = 0,l,2,.-- 

forftC-'^) eso(3,l), the Lie algebra of >S0(3,1). We assume X, = rX«+^~o r-'^X^-'^) 
and the time function of the embedding X^ is given by 



E 

fc=0 
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Therefore we may assume X*-^-* = (0, Xi,X2, is a standard embedding of S'^ into 
M^. We shall show that all f*^°\ ■ ■ ■ , f^"'^^ and b^^\ ■ ■ ■ , h^~^'> can be solved inductively. 
To prepare for the induction, we compute the corresponding terms in the optimal 
isometric embedding for of the given form fl4.8p . 
Fix an /, we denote 

X,+i = e'^-'^^-'^--e^-''''-\''""x. (4.9) 

Suppose and all h'^^\ ■ ■ - b^ are known, we see that other components 

of ■ ■ ■ are known and by (gSD all ■ ■ ■ are known. 

It is not hard to see that X^'^ = Xi~^ = ■ ■ ■ is stabilized at the r~' term and 
thus = Xi~2 . By definition, X1+2 is 

{/ + ■■■ + r-'-i[fo(-'-i) + G{b^'^\ ■ ■ ■ , M-'))] + ■ ■ ■ }(rBl(i) + ■ ■ ■ + r-'5l(-') + ■ ■ ■ ) 
= rSX« + ■ ■ ■ + r-'(b(-'-i)SX« + + G{b^'\ ■■■ , b^"'), f(°\ ■ ■ ■ , f(-'+i))) + 0{r~'-'), 

and thus 

and Xj^i^ for > m > -/ + 1 is of the form ■ ■ , M"'), fW, ■ ■ ■ , f(-'+i)). 

Denote by r the time function of Xi^2- We shall plug r into the optimal equation, 
find 6^"'"^) so that f*^"'^^'' is solvable, and then solve for f'^~'~^^ Denote by ©^'^^ a 
term of order r'^ whose coefficients depend only on f*^*^) , ■ ■ ■ , f'-"'"^^^ and b^^^ , ^'-"'^ 
and the physical data a, \H\, and V. 

Lemma 6. Write e'''°' = {3^^), then r = {Xi^2)o has the following expansion: 

r = rr^^^ + r-'r^"') + + ■ ■ ■ + (5^~^+^^ + 0{r-^~^) (4.10) 

where 

1=1 

and 

^^"^^ = E E C"'^^^^^^ + ^oof (-'^ + G(M°), . . . , M-'), f(°), ■ ■ ■ , f(-'+^)) (4.11) 
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Proof. We compute 
and recall (-^)o^'' = and {X)f'^ = Xi for i = 1,2,3. Likewise, 



On the other hand {BX^~'-^)q = Sqo'?"*-"''' + BQiX\~''\ and we already know that 
depends only on f^"'"'"^-*. □ 

Now we proceed to calculate the terms in the optimal isometric embedding equa- 
tion. For simplicity, we denote Y^i=iYl^a=Q^oa~^^ ^aiXi by g. 

Lemma 7. With r given by fl4.10p . we have 



At = -2r- V(i) - 2r~^~^g + (5^"^^ + ■■■ + ^^-'-'^ + 0{r~^-'^) 



and 



i i 

where 

a = Boi (4.12) 

3 

di = J2Ba^b'o-J-'^ (4.13) 

Proof. We use the formula that if two functions A and B on S"^ are given by A = AiXi 
and B = BjXj, then ■ Vfi = ^ AiBi - AB. □ 

Recalling the function / defined in (14.41) . we compute the expansion of / and the 
expansions of terms that appear on the optimal isometric embedding equation (14. 5 p 
in the following: 

Lemma 8. Suppose \Hr\ = 2r"^ + r^^/i^r^^ + (3^'^^ H h (S^-'-^) + ©(r"'""^) and 

\Ho\= 2r"i + r-2/i[,-2) + ©M) + . . . + 0(-«-3) ^ o{r^^~'^) then 

f = + + c5(-3) + . . . + + o(r-'-^) 
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where 

, (-2) _ , (-2) 

f (-2) = ^ ^ (4 14) 

J (l + |c|2)l/2 ^^-^^^ 



and 



^ (l + |cP)3/2 ^ ^ 1+|CP^" ^ ^ 



+ ©(-4) + . . . + 0(-«-4) ^ 



r '-5) 



4.3 Solving the optimal isometric embedding equation for 
all orders 

Theorem 2. Suppose satisfies (1), (2), and (3) and the limit of the quasilocal 
energy-momentum {e,pi,p2,P3) is timelike. There is a function t^'^^ on such that 
isometric embeddings : E — j- M^'^ with the time function r given below solves 



equation (JjJ^ W to 0{r 

3 

T = {Y,^,X,)r + T^''^ + 0{r-') 



where (01,02,03) satisfy 



i=l 



= = Si 

2 e 



Proof. Under the assumption, the energy momentum vector {e,Pi,P2,P3) and given 
by Equation (2.6)] 

/ (4~^^ - h^-^^)dS^ = Sttc and - / {vl~^^ - v^-^^)X^dS^ = Snp^. (4.16) 
From LemmalHl the coefficient of the leading r^~^^ term of div{fVT)—A sinh~^ \h^\\Ho\ 

is 

^(/(-2)vr(i))-lA(/(-2)ArW), 



while the '^^ term of divVr — divVo is 

^(-3)_lA(A + 2)f(°). 
2 

Thus f(°) is solvable if 

is zero. Computing using Vr^^^VXj = Q — T^^^^Xj, we find this is equivalent to 



52 



From fgrrel) we see that /^^ f-'^US'^ = -7=^8% and /^^ i-^-^^Xid^^ = 8np,. Thus 



we can take 

Ci Pi 



'l + |c|2 e 

and f(o) is solvable. By equation fl4.12p . b^'^^ is determined. □ 

Now we solve the optimal embedding equation to all order of r by induction. 

Theorem 3. Under the same assumption as in Theorem Suppose the family 
of optimal isometric embeddings X^. into 'E?''^ is of the form (14. 8 P described in the 
previous subsection. There exists b^~''^ G so(3, 1) such that all f^~^^ can be solved for 
k>0. 

Proof. From Lemma[8], the coefficient of the r'-"^"^-* term of div{fVT)—A sinh~^ \itl\Ho\ ' 
div{V — Vq) is 

rf^(/(-2) + y.H-3)y^(i)) _ 1 A(/(-'-3) A^(i) ^ fi~^)Ag) + i A(A + 2)f(-'-^) + e 

where is a term that depend on f^'^\ ■ ■ ■ f^~''\ b^^\ ■ ■ ■ b^~^\ a, V and H. 

Thus the solvability depends on whether, for each i = 1,2,3, we can make the 
following expression equal to — J^a dS"^ by choosing suitable dj in g: 



'52 



XidS^. 
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We integrate by parts and compute that this expression is equal to 



52 



where we use VgVXi = di — gXi 
This term is 



Since the energy component is positive and the matrix 5ij — is positive 

definite, we can choose dj so that Vz. 



(SX- dS' 



□ 



4.4 Locally energy minimizing at all orders 

In this subsection, we show the solution = B^X^ obtained in Theorem [3] is locally 
energy-minimizing at all orders. By this we mean that the second variation of the 
quasilocal energy is positive if we vary any f^"'^^ or 6^"^^^ Let r be the time function 
of X^ and let 5t be a variation. From Proposition 6.2 in [S], the first variation of the 
energy is 



-[11 a —a a 11, 



cd) 



divY.^{- 



Vr 



v/l + |Vr|2 v^l + |Vr|2 

By the derivation in §11 we can rewrite this as 

1 (A^)/ 



cosh^li/l -Ve-Vr 



drdJ^r 



s. 



div^fVr) — A sinh 



l-f^rll-f^ol 



dtv{V - Vo) 



drdHr- 



From Lemma [HI the coefficient of the r ' ^ term of div{f'VT) — A sinh ^ ^^^^-f 



div{V — Vo) is 



\Hr\\Ho\ 



+ ^A(A + 2)f(-'-i) + ■ • ■ f(°). 
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where g = J^Li ELo ^oJ B^t^i- 

On the other hand, the r~' term of r is 



r 



(-0=^ + 5QQfH) + G'(6W,-- - ,6(-'),f(o),--- ,f(-'+i)). (4.17) 



When we consider the variation of f*^ '\ 5f*^ '\ the leading term term of the 
second variation is of the order r~^'~'^ with coefficient 



00 



S2 



[A(A + 2)6f^~'\5f^~'^dS^. 



We may assume J^a Sf^"^^dS'^ = by normahzation. By decomposing 6f^~''^ into sum 
of eigenf unctions of S"^ and noting that the first non-zero eigenvalue of S"^ is —2, this 
is always positive. 

Varying is equivalent to varying g and 5g = ZlLi ELo '^^Oa'"^^^m-^i- 

The coefficient of the leading r~^'~^ term of the second variation with respect to 
^(-'-i) is then _ _ 

Using integration by part, we derive that this is equal to 

- / f^~'\\V6g\' + {Sgr]dS'. 
This simply gives a positive multiple of quasilocal energy since 

\V6g\' + {6gr = J2&i-J-'^B^^)" 

i a 

is a positive constant and — f^o f^~'^^dS'^ = Z^*^ . 
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